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ABSTRACT 


During the assembly flights of the Space Station Freedom (SSF), the Orbiter will ei- 
ther dock with the SSF and retract to the final berthed position, or will grapple the SSF 
using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled 
vehicles to their final berthed position. The SRMS method is expected to take approxi- 
mately one to one and a half hours to complete and require periodic attitude corrections 
by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a 
control moment gyro (CMC) system with RCS desaturation as required. Free drift of the 
attached vehicles is not currently thought to be acceptable because the desired system atti- 
tude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques 
resulting in power generation problems, thermodynamic control problems, and communi- 
cations problems. This paper deals with the simulation and control of the SRMS during 
tr un nion /latch interaction dynamics and during RCS maneuvers. The SRMS servo drive 
joints have highly non-linear elastic characteristics which tend to degrade sensitive control 
strategies. In addition the system natural frequencies are extremely low and depend on 
the drive joint deflections and SRMS geometric position. The lowest mean period of oscil- 
lation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the 
final berthed position is approximately 120 seconds. A detailed finite element model of the 
SRMS has been developed and used in a newly developed SRMS systems dynamics simu- 
lation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF 
systems. The present SRMS control strategy of brakes only recommended by the Charles 
Draper Labs is contrasted with a robust controller developed by the authors. The ro- 
bust controller uses an optimal linear quadratic regulator (LQR) to optimally place the 
closed-loop poles of a multivariable continuous-time system within the common region of 
an open sector with the sector angle ±45° from the negative real axis, and the left-hand 
side of a fine parallel to the imaginary axis in the complex s-plane. This guarantees that 
the critical damping ratio for the desired control modes is equal to or in excess of 0.707. 
The matrix sign function is used for solving the Riccati equations which appear in the 
controller design procedure. Fast and stable algorithms have recently been developed for 
the computation of the matrix sign function (ref. 1). Simulation results are given which 
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demonstrate the potential CSI involvement for the current SRMS control system and the 
proposed control system. 


INTRODUCTION 


The Space Station Freedom (SSF) will be assembled and supplied by the Orbiter. 
Each visit to the SSF will require a rigid structural attachment between the vehicles in 
order to assemble the mission build (MB) segments, and/or to resupply and exchange crew 
members. During the assembly flights of the SSF, the Orbiter will either dock to the SSF 
and retract to form a rigid and stable structural connection between the docking vehicles, 
or an Orbiter crew member will grapple the SSF using the SRMS and maneuver (see 
figure 1) to the berthed position. An intensive engineering study involving many Johnson 
Spacecraft Center (JSC) divisions, their subcontractors, together with the Canadian SPAR 
Aerospace Limited, is being conducted to determine the viability of the SRMS berthing 
approach for attaching the large space structures. The viability of the docking approach 

to mating vehicles with similar weight properties has been established in previous NASA 
projects. 

The SRMS attachment method for MB2 through MB6 is expected to take from one to 
one and a half hours to complete and requires the Orbiter to maneuver the SRMS attached 
system from the initial grapple attitude to the combined vehicle torque equilibrium attitude 
(TEA). Free drift of the attached vehicles at this point in the mission is not thought to 
be acceptable because the desired system TEA will quickly deteriorate to large attitude 
excursions due to the unbalanced gravity gradient and aerodynamic torques acting on the 
system. This will result in power supply problems, thermodynamic control problems, and 
communication problems. Alleviating these problems requires periodic attitude corrections 
by the Orbiter Reaction Control System (RCS), or the SSF RCS, or continuous control 
by the Control Moment Gyro (CMG) system with RCS desaturation as required. Where 
possible, the CMG system will be employed to maintain the TEA since its control is less 
likely to move the SRMS joints significantly. However, during the SRMS retraction while 
using the CMG’s it may be necessary to desaturate the moment build up of the CMG’s 
using either the Orbiter or the SSF RCS jets. The SRMS must not allow the berthing 
vehicles to contact each other with sufficient relative kinetic energy to damage any part of 
either structure. Thus for the SRMS attachment system to be viable, with regards to a 
structural loads standpoint, it must be able to withstand some form of attitude maneuvers 
using RCS jets, and be able to control and limit the relative contact kinetic energy and 
the constrained motion loading at the structural latch interfaces. 


The SRMS servo drive joints have highly nonlinear elastic characteristics due to the 
anti-backlash mechanisms in their gear boxes which tend to degrade sensitive control 
strategies. In addition the natural frequencies of the combined system are extremely low 
and nonlinear, depending on the drive joint deflections, the SRMS geometric position, 
and the control system gains if the brakes are not applied. The lowest mean period of 
oscillation for the 0rbiter/SRMS/MB6 system during the SSF trunnion insertion into the 
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Figure 1. Orbiter/SRMS/MB6 berthing orientation. 

Payload Retention Latch Assembly’s (PRLA’s), a relatively stiff SRMS position, is around 
120 seconds with the brakes on. The current SRMS system responds to external loads 
with non-commanded motion and long transient response oscillations which add critical- 
ly needed time to the retraction maneuver task. Evidence of this effect has been noted 
by crew members on past shuttle flights where crew members waited for the oscillation to 
decay before issuing new maneuver commands. The extremely low frequency response of 
the heavy Orbiter/SRMS/SSF system will only aggravate this aspect of the control task. 
A disturbance rejection methodology will greatly enhance SRMS operations. 

In order to accurately represent the structural dynamics of the entire coupled sys- 
tem and because of the complex nature of the SSF finite element component structures, 
a detailed finite element model of the SRMS with nonlinear joints was developed by the 
authors. This model has been coupled to the Orbiter and the finite element model of the 
MB6. The MB stage modal models are contractually furnished by the McDonnell Dou- 
glas Space System Company for the Space Station Freedom Loads and Dynamic Working 
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Group. Modal frequencies to 50 Hertz are retained for the analysis to accurately form a 
mathematical basis for the SRMS flexible body dynamics coupled with the joint servomo- 
tor dynamics. The requirements for loads convergence during trunnion to PRLA contact 
require modal frequencies to 25 Hertz. The coupled modal approach efficiently represents 
the arm segment to arm segment interaction since it is included in the basic finite 
element model with the correct boundary conditions. The modal model is then used in 
a newly developed SRMS systems dynamics simulation program, Manipulator Docking 
Dynamics Program (MDDP), to investigate the non-linear transient response dynamics of 
the Orbiter/SRMS/SSF systems during attitude control maneuvers and during interface 
trunnion/latch interaction dynamics. This approach has the advantage of accurately rep- 
resenting the entire SRMS flexural dynamics coupled with the flexible body dynamics of 
the attached vehicles. 

The present brakes on control system response of the SRMS has been simulated for a 
representative Orbiter PROS jet firing. A single .080 second pulse dual jet +pitch maneuver 
was executed and the system was allowed to ring out for 100 seconds. Next a similar case 
was simulated with the optimum robust controller developed by the authors. The robust 
controller uses a linear optimal quadratic regulator to optimally place the closed-loop poles 
of a multivariable continuous-time system within the common region of an open sector, 
and the left-hand side of a fine parallel to the imaginary axis in the complex s-plane. The 
open sector, see figure 2, is defined by a sector angle of ±45° from the negative real axis. 
This guarantees that the damping ratios for the desired control modes are in excess of 
0.707 percent of critical. The matrix sign function (ref. 1) is used for solving the Riccati 
equations which appear in the controller design procedure. 

In summary, this paper deals with (1) the development and capability of the analytical 
test bed MDDP used to simulate the systems dynamics and resulting transient structural 
dynamics, (2) the theoretical basis for an optimal quadratic regulator with pole placement 
controller developed by the authors and used as an optional control method in the SRMS, 
and (3) the simulation results of controlling the SRMSwith twodifferent controllers during 
the transient dynamics occurring when the PROS fires jets to begin a maneuver. 

SIMULATION DESCRIPTION 


A detailed finite element model of the SRMS has been developed to support the de- 
cision to study berthing the MB elements of SSF to the Orbiter using the SRMS. This 
model is intended to serve as a test bed in investigating the overall system dynamics dur- 
ing attitude control firing transients and during the interaction dynamics of the berthing 
vehicles as the SSF trunnions contact the Orbiter PRLA’s. The interaction contact dynam- 
ics have been correlated with result from a NASTRAN simulation. The NASTRAN run 
simulated an initial impact velocity between the left front PRLA and the corresponding 
SSF trunnion. The interface between the trunnion and PRLA was elastically represent- 
ed with a high stiffness spring perpendicular to the PRLA latch surface and with weak 
springs orthogonal to the load direction. Only the first compression loading segment of the 
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Figure 2. The region of interest in the continuous-time s-plane. 


NASTRAN response simulation is applicable for correlation of impact and bounce from 
the surface. The MDDP has the simulation capability to bounce from surface contact 
to surface contact while dynamically coupled with the SRMS system. Surface contact is 
maintained as long as there is a compressive load on the surface. Figures 3 and 4 give the 
resulting load history for the case of impacting the systems with a relative z 0 velocity of 
0.1 fps. In figure 3, the NASTRAN load is perpendicular to the surface, while in figure 4 
the MDDP load is given in component form. The bouncing characteristics predicted by 
MDDP are shown in figures 5 and 6 for a case where the trunnion first impacts near the 
apex of the latch then bounces from surface to surface. The load history is given m figure 
5 and the motion history is given in figure 6. The MDDP simulation will also be used to 
study the interaction dynamics of pre-integrated truss (PIT) to PIT berthing dynamics. 
The detailed transient load histories will be processed in load indicator subroutines to 
determine design limit loads for the SSF. The finite element model of the SRMS is cou- 
pled with a rigid body Orbiter and a finite element model of the SSF-MB6. The MDDP 
simulation of the SRMS includes the following capabilities: 

• Physical and dynamical properties 

• High frequency arm flexural dynamics (162 degrees of freedom ). 

• SRMS joint servo control model adapted from the JSC Shuttle Engineering Sim- 
ulator (SES), correlation in progress. 
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• Nonlinear hysteresis SRMS joints. 

• Second order servomotor dynamics. 

• SRMS brake dynamics including motor friction. 

• Locked or bottomed SRMS servo joints. 

• Flexible body payloads (716 degree of freedom for MB6 ). 

• Control Modes 

• Manual single joint drive with remaining joints in position hold. 

• Direct drive mode with remaining joints in brakes hold. 

• Position hold submode. 

• Optimal linear quadratic regulator with pole placement control system. 

The RCS attitude control systems for the Orbiter and the SSF are under development. 
The current program has the capability to apply a predetermined RCS pulse history to 
either vehicle for testing maneuver disturbances. 

The equations of motion for MDDP are linearized about the initial unstressed position 
geometry. For geometry changes greater than some e, the modes and frequencies of the 
system equations for the SRMS will be recomputed and the new modal state initialized 
such that the physical states of the total system match at the transition. For a fifty 
foot SRMS, a 5 percent geometry change would correspond to end deflections of 2.5 feet. 
Practically, modal recomputations do not occur very often for the type of studies that 
the MDDP is designed to handle. The relative motion interaction dynamics between the 
trunnions and PRLA’s are constrained by the latch geometry, thus limited to less than 
±3 inches in the Orbiter x 0 and y 0 directions. The Orbiter +z 0 direction is unlimited; 
however, the bottom of the PRLA latch limits the -z 0 motion direction. The top of the 
PRLA guides are a maximum of 24 inches tall. For the case of RCS excitement, a 2.5 foot 
change in geometry indicates that the system is not performing very well in its position 
hold capability. 

The flexible body model of the SSF was constructed by a component mode synthesis 
method using the Craig-Bampton (ref. 2) constrained boundary component mode reduc- 
tion technique. This dynamic reduction technique was used in both MSC/NASTRAN 
internally defined superelements and external superelements supplied by the different 
Space Station Work Package contractors and by the international partners. 
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The simulation is written in the FORTRAN 77 language and is computer based in 
the Engineering Computation Facility (ECF) at JSC. The ECF consists of a Cray X-MP 
EA/464 high speed vector processor (HVSP) and an Amdahl 5990-500 scalar processor, 
both operating under the UNIX operating environment. The simulation program is typi- 
cally submitted from the Amdahl front end to the Cray in batch form. 

In the following section, the optimal quadratic-pole placement technique used in this 
study will be presented together with an adaptation to the SRMS control system. 

OPTIMAL QUADRATIC REGULATOR WITH POLE PLACEMENT 


This section presents the mathematical basis for the optimal quadratic regulator as 
developed by the authors, and the adaptation to the SRMS system. The flexible body 
system equations of motion for the Orbiter/SRMS/SSF structure can be represented by 


Mx + Di + Kx = Cu; x(0), i(0) 


(1) 


where x and u are an n X 1 physical vector and an m X 1 control input vector, respectively, 
and M, D, K and C are of appropriate dimensions and represent the structural mass 
matrix, the damping matrix, the stiffness matrix, and the control distribution matrix, 
respectively. 

Next we employ modal analysis to Eq. (1) to decouple and reduce the order of 
the problem to a smaller size based on the dominate modes. The physical to modal 
transformation is given by 


x = 4 > y, 


( 2 ) 


where <f> is the n X n mode shape matrix and y is the n X 1 modal vector. Using Eq. (2) 
in Eq. (1) yields the decoupled modal equation (ref. 3) given by 

Iy + 2£u>y + u> 2 y = ^ T Cu. (3) 

where I = ^^M^; and ti> 2 = The superscript T denotes the 

transpose of a matrix. 


We reorder our modal matrices to retain the elements we wish to design in the optimal 
control section in the following manner 
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(4) 


y = 


z r 

Zo 


where z r are the regulated modal variables to be designed and z 0 are the omitted modal 
variables. We then obtain the following 


Z+ 
Z o 


‘2^ 7 4V r 

0 


0 

2 ^ 0^0 






(5) 


where z r is an n r X 1 vector, with appropriate dimensions for the remaining terms. The 
equations to design by optimal controller can be separated from Eq. (5) as 


z,. -j- "f n> r z r — C,u, 


( 6 ) 


Next we can express Eq. (6) in the first order state space form in terms of a new 
vector w as 


Wi 


w 2 



0 


it} 


I 

2(rU} r 


Wi 

W 2 


+ 


0 

CL 


u. 


(7) 


Equation (7) in compact notation becomes the linear controllable continuous-time 
system described by 


w = Aw + Bu, 


w( 0 ) 


(8) 


where w and u are 2 n r x 1 and m x 1 respectively. A and B follow from Eq. (7) and 
are constant matrices of appropriate dimensions. Let the quadratic cost function for the 
system in Eq. (8) be 


>• + 00 

J = / [w t Qw + u T Ru]dt, ( 9 ) 

where the weighting matrices Q and R are 2n r x 2n r non-negative definite and m x m 
positive definite symmetric matrices, respectively. The feedback control law that minimizes 
the performance index J in Eq. (9) is given by (ref. 4) : 
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u = -Kw + Er = — R -1 B t Pw + Er 


( 10 ) 


where K is an m X 2 n r feedback gain; E is an m x m forward gain; r is an m X 1 reference 
input; and P, an 2n r X 2 n T non-negative definite symmetric matrix, is the solution of the 
Riccati equation, 


PBR -1 B t P - PA - A t P - Q = 0 (11) 

with (A, Q) detectable. Thus the resulting closed-loop system becomes 


w = (A — BK)w + BEr 


( 12 ) 


The eigenvalues of (A — BK), denoted by cr(A — BK), he in the open left— half plane 
of the complex s plane. Our objective given a suitable R is to determine Q, P, and K so 
that the closed-loop system in Eq. (12) has its eigenvalues on or within a specified sector 
region shown in figure 2. The important results to achieve the desired design are presented 
in the following. 

Lemma 1. : (refs. 4 and 5) Let (A, B) be the pair of the given open-loop system in Eq. 
(8) with the quadratic cost function in Eq. (9). Also let h > 0 represent the prescribed 
degrees of relative stability. Then the eigenvalues of the closed-loop system (A — BK) 
consist of the invariant eigenvalues of A lying to the left of the —h vertical line with the 
matrix P being the solution of the Riccati equation: 


PBR _1 B r P - P(A + hi) - (A + hI) T P - 0 (13) 

Theorem 1 : (ref. 6) Let the given stable system matrix A 6 5i nX " have eigenvalues 
cr7(i — 1,2, .. . ,n~) lying in the open sector of figure 2 with the sector angle ±45° from 
the negative real axis and the eigenvalues (z = 1,2, . . . outside that sector, with 

n = n“ + n + . Now, consider the two Riccati equations, 

QBR _1 B r Q - Q(-A 2 ) - (A 2 ) t Q = 0 (14a) 


and 


PBR -1 B t P - PA - A t P - Q = 0 


(146) 
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Then, the closed-loop system, 


A c = A- 7 BK = A- 7 BR- 1 B t P, (15) 

will enclose the invariant eigenvalues a~(i = 1 and at least one additional pair 

of complex conjugate eigenvalues lying in the open sector of figure 2, for the constant gain 
7 in Eq. (15) satisfying 


7 > max{ 


1 b + ^/b 2 + ac 

, / 
z a 


( 16 ) 


where . = lr[(BR-'B T P) ! |, b = (r[BR-‘B T PA], and c = (l/2)(r[BR->B T Q|. 

Remark 1 : The steady state solutions of the Riccati equations in Eqs. (11 and 14) can 
be found using the matrix sign function techniques (refs. 7, 8). 

The steps to optimally place all the closed-loop eigenvalues in the hatched region of 
figure 2 are described as follows. 

Step 1 : Let the given continuous time system be as in Eq. (8). Specify h so that the — h 
vertical line on the negative real axis would represent the line beyond which the eigenvalues 
have to be placed in the cross hatched sector of figure 2. Also, assign Ao = A and the 
positive definite matrix R = I. Set i = 1. If the system is unstable, then solve Eq. (13) 
to obtain P 0 and the closed-loop system Ai — A — 7 0 BR~ 1 B T P 0 = A — 7 0 BK 0 , with 
7 o = 1; else (stable system) go to Step 2 with Ai = A,P 0 = 0 and 70 = 0 . 

Step 2 : Solve Eq. (14a) for Q, with A := A,-. Check if (l/2)tr[BR _1 B T Q t -] is zero. If 
it is equal to zero, go to Step 4 with j=i; else, continue and go to Step 3. Note that, when 
(l/2)fr[BR _1 B T Q,] = 0, all eigenvalues of the matrix A,* He on or within the open sector 
of figure 2. 

Step 3 : Solve Eq. (14b) for P, with A := A; and Q Q t -. Then the constant gain 7 ,- 
can be evaluated using Eq. (16). The closed-loop system matrix is 


A 1+] = Ai - 7,BR 1 B t P 1 - = Ai - 7 ,BK, (17a) 

Set i i + 1 and go to Step 2. 

Step 4 : Check if <r[(A j + fol)] + (sum of the eigenvalues to the right of the vertical line at 
— h) is zero. If it is equal to zero, go to Step 5 with Pj+i = 0 and 7 j + i = 0; else, solve Eq. 
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(13) for P j+ i with A := A j and obtain the closed-loop system Aj -7j+iBR 1 B T Pj + i - 
A j - 7 j+iBKj+i, with 7 j+i = 1 and K j+1 = R _1 B T Pj+i. 

Step 5 : The designed closed-loop system is 


j + 1 

A 0 - BR 1 B T J2^ P k ( 176 ) 

k = 0 

and its eigenvalues he in the hatched region of figure 2. Note that the above system matrix 
in Eq. (17b) is equal to the system matrix in Eq. (12), with 


i 

Q = 2fcP 0 + £(Qi + A 7l P i BR- 1 B T P l ) 7 ,- (18) 

1=1 


In the above equation, A 7 , - = 7t - — 1 and the matrix R is as originally assigned. Also, 
the optimal continuous time regulator can be given as 

j + 1 

)w + Er = — Kw + Er (19) 

i= o 

where r is any reference input and K is the desired state feedback gain matrix. 

Next we transform back to the modal state y in order to apply our control gains back 
in our test bed simulation. Note that control input u becomes 


u = — Kw = — [Kn K 


12 


Z •p 

z r 


( 20 ) 


and further rearrangement yields 


u = — [Kn 0] 


z r 

z 0 


-[K 12 0] 


Z r 

io 


( 21 ) 


and using Eq. (4) we obtain the desired input 


u = — [Kn 0]y — [Ki 2 0]y. 


( 22 ) 
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In the test bed we solve for u and substitute back into Eq. (3) and use a variable step 
fourth order Runge Kutta integrating method to propagate a time history of the system 
dynamics. 

The control system feedback gain matrix, K, was obtained by following the above de- 
sign steps for the case of the Orbiter SRMS handling the MB6 spacecraft. The complex fre- 
quencies for the dominate modes of the open loop system with soft torsional springs at each 
of the SRMS servo joints are { -.000153 ± >.00766; -.000235T >.0112; -.000595 ± >. 0297; 

— .000923 ± >.0461; —.00128 it >.0638; —.00535 ± >.268 }. The value of h = .004 was 
chosen to place the closed loop poles in the phase plane such that the joint torque require- 
ments u do not exceed the system capability for the given PRCS disturbance. The closed 
loop system eigenvalues tr(A-BK) after design are { -.00785 ± >.00766; -.0145 ± >.118; 

— .054 ± >.0468; — .103 ± >.0937; — .166± >.157; — 1.594 ± >1.594 } all within the required 
stability sector given in figure 2. The feedback gain matrices for this design are given as 


and 


[Kn] = 


N = 


■ — .2981el 
+.1588el 

— .2854e0 
+ .1249e2 

— .2795el 
. +.5814el 


-+.3766e3 
+.3493e3 
— .3622e3 
+.2473e4 
+.7537e3 
. — .1587e4 


+.2026el 
+ .1997e0 

— .5697e0 

— .4620e0 
-f .2346el 

— .5165el 


-f .2269el 

— .1078e3 
+.4103e2 

— .6692e3 
+.7676e2 

— .1212e3 


— .5281e0 

— .2833el 
+ .1168e2 

— .2918el 
+ .1060e2 
+.6126el 


— .4313e2 

— .1035e2 
+ .1935e3 
+.2952e2 
+.3337e3 
+ .1609e3 


+.3582el 

— .3350e2 
+.1785e2 
+.6502el 

— .5890el 

— .8879el 


+.1937e2 

— .3947e3 
+.1719e3 
+.1214e3 
— .3160e2 

— .9004e2 


+.2090e2 

— .8016el 
+.5723e2 

— .5299e2 

— .4631e2 
+.1820e2 


+.1539e3 

— .4719e2 
T .4239e3 

— .3278e3 

— .2773e3 
+.1312e3 


-,8834e3- 

— .1452e2 
+ .1112e2 

— .1330e4 
+.4805e3 

— .1574e4. 


-,5734e3- 

— .1115e2 
+ .1438e2 

— .8417e3 
+.3002e3 

— .9929e3. 


SIMULATION RESULTS 


The SRMS was designed to handle payload weights up to 45,000 lbs. with the Orbiter 
attitude control system in a free drift mode. Current plans to grapple vehicles such as SSF- 
MB6, which weighs approximately 250,000 lbs., and to maneuver the combined spacecraft 
require a close engineering look at the systems and careful testing. The new alternative 
digital auto pilot (ALT DAP ) control mode for the Orbiter, which wasn’t available when 
the SRMS was originally designed, will greatly help reduce the dynamic excitation of the 
system. The ALT DAP allows both selectable burn times and delay times between firings, 
so conceivably the elastic response disturbances can have a chance to damp out between 
firings. With this in mind, a test case representing a .08 second Orbiter PRCS dual jet 


522 


pitch axis firing was simulated to determine the dynamic response motion between the 
Orbiter PRLA’s and the MB6 trunnions. 

The first case was simulated using anSRMS brakes hold mode for all servo joints with 
the arm geometry placing the pressurized berthing adapter (PBA) in the cargo bay of 
the Orbiter with the SSF trunnions inside the PRLA’s. The interaction loads dynamics 
between the trunnions and the PRLA’s was turned off in MDDP m order to see the 
unrestrained motion. The case was simulated for 100 seconds with the initial state of 
the system being zero. Shown in figures 7, 8, and 9 are the Orbiter’s angular rate vector 
history together with the total system angular rate vector history, the velocity history of 
the left front trunnion with respect to the left front PRLA, and the motion of the left 
front trunnion with respect to the left front PRLA in a x 0 ,— y 0 plane, respectively. Figure 7 
shows that the system damps very slow after the brakes quit slipping. This represents a poten- 
tial problem, since additional thrusting can easily reinforce the excitation level, causing 
additional slip of the motors and payloads. The relative velocities shown in figure 8 are less 
than .013 fps after 20 seconds. The relative motion shown in figure 9 did not exceed 

The second case simulated used the optimal controller for 100 seconds of simulation 
time. The results shown in figures 10, 11, and 12 are for the same parameters as in the first 
case. The angular rate results shown in figure 10 exhibit the optimal characteristic of a 
slight dynamic overshoot of approximately 0.008 dps for the Orbiter roll axis, and remain 
below this value after about 17 seconds as the controller brings the Orbiter rates to match 
with the overall system rates. The relative velocities between the trunnion and PRLA 
in figure 11 stay below 0.005 fps after 12 seconds from the initial transient PROS firing. 
Very small amplitude modal coupling can be seen in the results for the optimal controller 
as compared to the coupling seen in previous cases. The optimal controller quickly rejects 
the PROS disturbance, reducing the system state to a level where additional firings can 
be made. The motion of the trunnion with respect to the Orbiter shown in figure 12 was 
very smooth and peaked at 0.26 feet, slightly more than for the case with applied brakes. 

CONCLUSIONS 


This paper has presented a description of MDDP, a high fidelity finite element model 
simulation of the SRMS specifically designed to study SRMS structural transient responses 
while connected to a flexible body payload. Of particular interest to structural engineers 
are the interaction dynamics computed between structural elements in the simulation and 
also the familiar finite element form of the equations. The simulation addresses two specific 
dynamics events, (1) constrained motion transient dynamics and the resulting structural 
loads occurring during interface contact, and (2) system transient dynamics resulting from 
the firing of RCS jets. The constrained motion loading relies entirely on the finite element 
representation of the contacting structures since no attenuation mechanism has been allo- 
cated for this interface. A new optimal LQR control method has been presented and shown 
to be very effective for rejecting the transient responses caused by the firing of the Orbiter 
PROS. Future work will test the optimal LQR controller during the interaction dynamics 
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contact phase. The multivariable state- space optimal controller is stable and has good 
robustness properties. The state-space approach mathematically decouples the actuators 
from each other thereby greatly reducing the unwanted dynamic motion typically seen in 
the SRMS transient response. 
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Figure 9. Trunnion motion in the xo-zo plane with respect to the PRLA for brake hold. 
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Figure 10. Orbiter and Orbiter/SRMS/MB6 angular rate history for optimal control. 
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Figure 12. Trunnion motion in the xo-zo plane with respect to the PRLA for optimal control 









